Multimodal machine learning for modeling infant head circumference, mothers' milk composition, and their shared environment.

Links between human milk (HM) and infant development are poorly understood and often focus on individual HM components. Here we apply multi-modal predictive machine learning to study HM and head circumference (a proxy for brain development) among 1022 mother-infant dyads of the CHILD Cohort. We integrated HM data (19 oligosaccharides, 28 fatty acids, 3 hormones, 28 chemokines) with maternal and infant demographic, health, dietary and home environment data. Head circumference was significantly predictable at 3 and 12 months. Two of the most associated features were HM n3-polyunsaturated fatty acid C22:6n3 (docosahexaenoic acid, DHA; p = 9.6e-05) and maternal intake of fish (p = 4.1e-03), a key dietary source of DHA with established relationships to brain function. Thus, using a systems biology approach, we identified meaningful relationships between HM and brain development, which validates our statistical approach, gives credence to the novel associations we observed, and sets the foundation for further research with additional cohorts and HM analytes.

Large-scale correlation network construction for unraveling the coordination of complex biological systems.

Advanced measurement and data storage technologies have enabled high-dimensional profiling of complex biological systems. For this, modern multiomics studies regularly produce datasets with hundreds of thousands of measurements per sample, enabling a new era of precision medicine. Correlation analysis is an important first step to gain deeper insights into the coordination and underlying processes of such complex systems. However, the construction of large correlation networks in modern high-dimensional datasets remains a major computational challenge owing to rapidly growing runtime and memory requirements. Here we address this challenge by introducing CorALS (Correlation Analysis of Large-scale (biological) Systems), an open-source framework for the construction and analysis of large-scale parametric as well as non-parametric correlation networks for high-dimensional biological data. It features off-the-shelf algorithms suitable for both personal and high-performance computers, enabling workflows and downstream analysis approaches. We illustrate the broad scope and potential of CorALS by exploring perspectives on complex biological processes in large-scale multiomics and single-cell studies.

Deficiency of the Heterogeneous Nuclear Ribonucleoprotein U locus leads to delayed hindbrain neurogenesis.

Genetic variants affecting Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU) have been identified in several neurodevelopmental disorders (NDDs). HNRNPU is widely expressed in the human brain and shows the highest postnatal expression in the cerebellum. Recent studies have investigated the role of HNRNPU in cerebral cortical development, but the effects of HNRNPU deficiency on cerebellar development remain unknown. Here, we describe the molecular and cellular outcomes of HNRNPU locus deficiency during in vitro neural differentiation of patient-derived and isogenic neuroepithelial stem cells with a hindbrain profile. We demonstrate that HNRNPU deficiency leads to chromatin remodeling of A/B compartments, and transcriptional rewiring, partly by impacting exon inclusion during mRNA processing. Genomic regions affected by the chromatin restructuring and host genes of exon usage differences show a strong enrichment for genes implicated in epilepsies, intellectual disability, and autism. Lastly, we show that at the cellular level HNRNPU downregulation leads to an increased fraction of neural progenitors in the maturing neuronal population. We conclude that the HNRNPU locus is involved in delayed commitment of neural progenitors to differentiate in cell types with hindbrain profile.

Deep representation learning identifies associations between physical activity and sleep patterns during pregnancy and prematurity.

Preterm birth (PTB) is the leading cause of infant mortality globally. Research has focused on developing predictive models for PTB without prioritizing cost-effective interventions. Physical activity and sleep present unique opportunities for interventions in low- and middle-income populations (LMICs). However, objective measurement of physical activity and sleep remains challenging and self-reported metrics suffer from low-resolution and accuracy. In this study, we use physical activity data collected using a wearable device comprising over 181,944 h of data across N = 1083 patients. Using a new state-of-the art deep learning time-series classification architecture, we develop a 'clock' of healthy dynamics during pregnancy by using gestational age (GA) as a surrogate for progression of pregnancy. We also develop novel interpretability algorithms that integrate unsupervised clustering, model error analysis, feature attribution, and automated actigraphy analysis, allowing for model interpretation with respect to sleep, activity, and clinical variables. Our model performs significantly better than 7 other machine learning and AI methods for modeling the progression of pregnancy. We found that deviations from a normal 'clock' of physical activity and sleep changes during pregnancy are strongly associated with pregnancy outcomes. When our model underestimates GA, there are 0.52 fewer preterm births than expected (P = 1.01e - 67, permutation test) and when our model overestimates GA, there are 1.44 times (P = 2.82e - 39, permutation test) more preterm births than expected. Model error is negatively correlated with interdaily stability (P = 0.043, Spearman's), indicating that our model assigns a more advanced GA when an individual's daily rhythms are less precise. Supporting this, our model attributes higher importance to sleep periods in predicting higher-than-actual GA, relative to lower-than-actual GA (P = 1.01e - 21, Mann-Whitney U). Combining prediction and interpretability allows us to signal when activity behaviors alter the likelihood of preterm birth and advocates for the development of clinical decision support through passive monitoring and exercise habit and sleep recommendations, which can be easily implemented in LMICs.

Cross-species comparative analysis of single presynapses.

Comparing brain structure across species and regions enables key functional insights. Leveraging publicly available data from a novel mass cytometry-based method, synaptometry by time of flight (SynTOF), we applied an unsupervised machine learning approach to conduct a comparative study of presynapse molecular abundance across three species and three brain regions. We used neural networks and their attractive properties to model complex relationships among high dimensional data to develop a unified, unsupervised framework for comparing the profile of more than 4.5 million single presynapses among normal human, macaque, and mouse samples. An extensive validation showed the feasibility of performing cross-species comparison using SynTOF profiling. Integrative analysis of the abundance of 20 presynaptic proteins revealed near-complete separation between primates and mice involving synaptic pruning, cellular energy, lipid metabolism, and neurotransmission. In addition, our analysis revealed a strong overlap between the presynaptic composition of human and macaque in the cerebral cortex and neostriatum. Our unique approach illuminates species- and region-specific variation in presynapse molecular composition.

Whole genome deconvolution unveils Alzheimer's resilient epigenetic signature.

Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) accurately depicts the chromatin regulatory state and altered mechanisms guiding gene expression in disease. However, bulk sequencing entangles information from different cell types and obscures cellular heterogeneity. To address this, we developed Cellformer, a deep learning method that deconvolutes bulk ATAC-seq into cell type-specific expression across the whole genome. Cellformer enables cost-effective cell type-specific open chromatin profiling in large cohorts. Applied to 191 bulk samples from 3 brain regions, Cellformer identifies cell type-specific gene regulatory mechanisms involved in resilience to Alzheimer's disease, an uncommon group of cognitively healthy individuals that harbor a high pathological load of Alzheimer's disease. Cell type-resolved chromatin profiling unveils cell type-specific pathways and nominates potential epigenetic mediators underlying resilience that may illuminate therapeutic opportunities to limit the cognitive impact of the disease. Cellformer is freely available to facilitate future investigations using high-throughput bulk ATAC-seq data.

Combined proteomics and CRISPR‒Cas9 screens in PDX identify ADAM10 as essential for leukemia in vivo.

BACKGROUND: Acute leukemias represent deadly malignancies that require better treatment. As a challenge, treatment is counteracted by a microenvironment protecting dormant leukemia stem cells. METHODS: To identify responsible surface proteins, we performed deep proteome profiling on minute numbers of dormant patient-derived xenograft (PDX) leukemia stem cells isolated from mice. Candidates were functionally screened by establishing a comprehensive CRISPR‒Cas9 pipeline in PDX models in vivo. RESULTS: A disintegrin and metalloproteinase domain-containing protein 10 (ADAM10) was identified as an essential vulnerability required for the survival and growth of different types of acute leukemias in vivo, and reconstitution assays in PDX models confirmed the relevance of its sheddase activity. Of translational importance, molecular or pharmacological targeting of ADAM10 reduced PDX leukemia burden, cell homing to the murine bone marrow and stem cell frequency, and increased leukemia response to conventional chemotherapy in vivo. CONCLUSIONS: These findings identify ADAM10 as an attractive therapeutic target for the future treatment of acute leukemias.

Screening autism-associated environmental factors in differentiating human neural progenitors with fractional factorial design-based transcriptomics.

Research continues to identify genetic variation, environmental exposures, and their mixtures underlying different diseases and conditions. There is a need for screening methods to understand the molecular outcomes of such factors. Here, we investigate a highly efficient and multiplexable, fractional factorial experimental design (FFED) to study six environmental factors (lead, valproic acid, bisphenol A, ethanol, fluoxetine hydrochloride and zinc deficiency) and four human induced pluripotent stem cell line derived differentiating human neural progenitors. We showcase the FFED coupled with RNA-sequencing to identify the effects of low-grade exposures to these environmental factors and analyse the results in the context of autism spectrum disorder (ASD). We performed this after 5-day exposures on differentiating human neural progenitors accompanied by a layered analytical approach and detected several convergent and divergent, gene and pathway level responses. We revealed significant upregulation of pathways related to synaptic function and lipid metabolism following lead and fluoxetine exposure, respectively. Moreover, fluoxetine exposure elevated several fatty acids when validated using mass spectrometry-based metabolomics. Our study demonstrates that the FFED can be used for multiplexed transcriptomic analyses to detect relevant pathway-level changes in human neural development caused by low-grade environmental risk factors. Future studies will require multiple cell lines with different genetic backgrounds for characterising the effects of environmental exposures in ASD.

Endothermic physiology of extinct megatooth sharks.

The evolution of the extinct megatooth shark, Otodus megalodon, and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark's gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

Evidence for human milk as a biological system and recommendations for study design-a report from "Breastmilk Ecology: Genesis of Infant Nutrition (BEGIN)" Working Group 4.

Human milk contains all of the essential nutrients required by the infant within a complex matrix that enhances the bioavailability of many of those nutrients. In addition, human milk is a source of bioactive components, living cells and microbes that facilitate the transition to life outside the womb. Our ability to fully appreciate the importance of this matrix relies on the recognition of short- and long-term health benefits and, as highlighted in previous sections of this supplement, its ecology (i.e., interactions among the lactating parent and breastfed infant as well as within the context of the human milk matrix itself). Designing and interpreting studies to address this complexity depends on the availability of new tools and technologies that account for such complexity. Past efforts have often compared human milk to infant formula, which has provided some insight into the bioactivity of human milk, as a whole, or of individual milk components supplemented with formula. However, this experimental approach cannot capture the contributions of the individual components to the human milk ecology, the interaction between these components within the human milk matrix, or the significance of the matrix itself to enhance human milk bioactivity on outcomes of interest. This paper presents approaches to explore human milk as a biological system and the functional implications of that system and its components. Specifically, we discuss study design and data collection considerations and how emerging analytical technologies, bioinformatics, and systems biology approaches could be applied to advance our understanding of this critical aspect of human biology.

Multiomic signals associated with maternal epidemiological factors contributing to preterm birth in low- and middle-income countries.

Preterm birth (PTB) is the leading cause of death in children under five, yet comprehensive studies are hindered by its multiple complex etiologies. Epidemiological associations between PTB and maternal characteristics have been previously described. This work used multiomic profiling and multivariate modeling to investigate the biological signatures of these characteristics. Maternal covariates were collected during pregnancy from 13,841 pregnant women across five sites. Plasma samples from 231 participants were analyzed to generate proteomic, metabolomic, and lipidomic datasets. Machine learning models showed robust performance for the prediction of PTB (AUROC = 0.70), time-to-delivery (r = 0.65), maternal age (r = 0.59), gravidity (r = 0.56), and BMI (r = 0.81). Time-to-delivery biological correlates included fetal-associated proteins (e.g., ALPP, AFP, and PGF) and immune proteins (e.g., PD-L1, CCL28, and LIFR). Maternal age negatively correlated with collagen COL9A1, gravidity with endothelial NOS and inflammatory chemokine CXCL13, and BMI with leptin and structural protein FABP4. These results provide an integrated view of epidemiological factors associated with PTB and identify biological signatures of clinical covariates affecting this disease.

Ribonucleoprotein (RNP)-Mediated Targeted Mutagenesis in Barley (Hordeum vulgare L.).

The crop species barley is a genetic model for the small grain temperate cereals. Thanks to the availability of whole genome sequence and the development of customizable endonucleases, site-directed genome modification has recently revolutionized genetic engineering. Several platforms have been established in plants, with the most flexible one offered by the clustered regularly interspaced short palindromic repeats (CRISPR) technology. In this protocol, commercially available synthetic guide RNAs (gRNAs), Cas enzymes, or custom-generated reagents are used for targeted mutagenesis in barley. The protocol has been successfully used with immature embryo explants to generate site-specific mutations in regenerants. As the double-strand break-inducing reagents are customizable and can be efficiently delivered, pre-assembled ribonucleoprotein (RNP) complexes allow efficient generation of genome-modified plants.

Ribonucleoprotein (RNP)-Mediated Allele Replacement in Barley (Hordeum vulgare L.) Leaves.

Varietal differences within a species with agronomic importance are often based on minor changes in the genomic sequence. For example, fungus-resistant and fungus-susceptible wheat varieties may vary in only one amino acid. The situation is similar with the reporter genes Gfp and Yfp where two base pairs cause a shift in the emission spectrum from green to yellow. Methods of targeted double-strand break induction now allow this exchange precisely with the simultaneous transfer of the desired repair template. However, these changes rarely lead to a selective advantage that can be used in generating such mutant plants. The protocol presented here allows a corresponding allele replacement at the cellular level using ribonucleoprotein complexes in combination with an appropriate repair template. The efficiencies achieved are comparable to other methods with direct DNA transfer or integration of the corresponding building blocks in the host genome. They are in the range of 35 percent, considering one allele in a diploid organism as barley and using Cas9 RNP complexes.

Data-driven longitudinal characterization of neonatal health and morbidity.

Although prematurity is the single largest cause of death in children under 5 years of age, the current definition of prematurity, based on gestational age, lacks the precision needed for guiding care decisions. Here, we propose a longitudinal risk assessment for adverse neonatal outcomes in newborns based on a deep learning model that uses electronic health records (EHRs) to predict a wide range of outcomes over a period starting shortly before conception and ending months after birth. By linking the EHRs of the Lucile Packard Children's Hospital and the Stanford Healthcare Adult Hospital, we developed a cohort of 22,104 mother-newborn dyads delivered between 2014 and 2018. Maternal and newborn EHRs were extracted and used to train a multi-input multitask deep learning model, featuring a long short-term memory neural network, to predict 24 different neonatal outcomes. An additional cohort of 10,250 mother-newborn dyads delivered at the same Stanford Hospitals from 2019 to September 2020 was used to validate the model. Areas under the receiver operating characteristic curve at delivery exceeded 0.9 for 10 of the 24 neonatal outcomes considered and were between 0.8 and 0.9 for 7 additional outcomes. Moreover, comprehensive association analysis identified multiple known associations between various maternal and neonatal features and specific neonatal outcomes. This study used linked EHRs from more than 30,000 mother-newborn dyads and would serve as a resource for the investigation and prediction of neonatal outcomes. An interactive website is available for independent investigators to leverage this unique dataset: https://maternal-child-health-associations.shinyapps.io/shiny_app/.

Prediction of neuropathologic lesions from clinical data.

INTRODUCTION: Post-mortem analysis provides definitive diagnoses of neurodegenerative diseases; however, only a few can be diagnosed during life. METHODS: This study employed statistical tools and machine learning to predict 17 neuropathologic lesions from a cohort of 6518 individuals using 381 clinical features (Table S1). The multisite data allowed validation of the model's robustness by splitting train/test sets by clinical sites. A similar study was performed for predicting Alzheimer's disease (AD) neuropathologic change without specific comorbidities. RESULTS: Prediction results show high performance for certain lesions that match or exceed that of research annotation. Neurodegenerative comorbidities in addition to AD neuropathologic change resulted in compounded, but disproportionate, effects across cognitive domains as the comorbidity number increased. DISCUSSION: Certain clinical features could be strongly associated with multiple neurodegenerative diseases, others were lesion-specific, and some were divergent between lesions. Our approach could benefit clinical research, and genetic and biomarker research by enriching cohorts for desired lesions.

Objective Activity Parameters Track Patient-specific Physical Recovery Trajectories After Surgery and Link With Individual Preoperative Immune States.

OBJECTIVE: The longitudinal assessment of physical function with high temporal resolution at a scalable and objective level in patients recovering from surgery is highly desirable to understand the biological and clinical factors that drive the clinical outcome. However, physical recovery from surgery itself remains poorly defined and the utility of wearable technologies to study recovery after surgery has not been established. BACKGROUND: Prolonged postoperative recovery is often associated with long-lasting impairment of physical, mental, and social functions. Although phenotypical and clinical patient characteristics account for some variation of individual recovery trajectories, biological differences likely play a major role. Specifically, patient-specific immune states have been linked to prolonged physical impairment after surgery. However, current methods of quantifying physical recovery lack patient specificity and objectivity. METHODS: Here, a combined high-fidelity accelerometry and state-of-the-art deep immune profiling approach was studied in patients undergoing major joint replacement surgery. The aim was to determine whether objective physical parameters derived from accelerometry data can accurately track patient-specific physical recovery profiles (suggestive of a 'clock of postoperative recovery'), compare the performance of derived parameters with benchmark metrics including step count, and link individual recovery profiles with patients' preoperative immune state. RESULTS: The results of our models indicate that patient-specific temporal patterns of physical function can be derived with a precision superior to benchmark metrics. Notably, 6 distinct domains of physical function and sleep are identified to represent the objective temporal patterns: ''activity capacity'' and ''moderate and overall activity (declined immediately after surgery); ''sleep disruption and sedentary activity (increased after surgery); ''overall sleep'', ''sleep onset'', and ''light activity'' (no clear changes were observed after surgery). These patterns can be linked to individual patients preopera-tive immune state using cross-validated canonical-correlation analysis. Importantly, the pSTAT3 signal activity in monocytic myeloid-derived suppressor cells predicted a slower recovery. CONCLUSIONS: Accelerometry-based recovery trajectories are scalable and objective outcomes to study patient-specific factors that drive physical recovery.

Deleterious and Protective Psychosocial and Stress-Related Factors Predict Risk of Spontaneous Preterm Birth.

OBJECTIVES: The aim of the study was to: (1) Identify (early in pregnancy) psychosocial and stress-related factors that predict risk of spontaneous preterm birth (PTB, gestational age <37 weeks); (2) Investigate whether "protective" factors (e.g., happiness/social support) decrease risk; (3) Use the Dhabhar Quick-Assessment Questionnaire for Stress and Psychosocial Factors (DQAQ-SPF) to rapidly quantify harmful or protective factors that predict increased or decreased risk respectively, of PTB. STUDY DESIGN: This is a prospective cohort study. Relative risk (RR) analyses investigated association between individual factors and PTB. Machine learning-based interdependency analysis (IDPA) identified factor clusters, strength, and direction of association with PTB. A nonlinear model based on support vector machines was built for predicting PTB and identifying factors that most strongly predicted PTB. RESULTS: Higher levels of deleterious factors were associated with increased RR for PTB: General anxiety (RR = 8.9; 95% confidence interval [CI] = 2.0,39.6), pain (RR = 5.7; CI = 1.7,17.0); tiredness/fatigue (RR = 3.7; CI = 1.09,13.5); perceived risk of birth complications (RR = 4; CI = 1.6,10.01); self-rated health current (RR = 2.6; CI = 1.0,6.7) and previous 3 years (RR = 2.9; CI = 1.1,7.7); and divorce (RR = 2.9; CI = 1.1,7.8). Lower levels of protective factors were also associated with increased RR for PTB: low happiness (RR = 9.1; CI = 1.25,71.5); low support from parents/siblings (RR = 3.5; CI = 0.9,12.9), and father-of-baby (RR = 3; CI = 1.1,9.9). These factors were also components of the clusters identified by the IDPA: perceived risk of birth complications (p < 0.05 after FDR correction), and general anxiety, happiness, tiredness/fatigue, self-rated health, social support, pain, and sleep (p < 0.05 without FDR correction). Supervised analysis of all factors, subject to cross-validation, produced a model highly predictive of PTB (AUROC or area under the receiver operating characteristic = 0.73). Model reduction through forward selection revealed that even a small set of factors (including those identified by RR and IDPA) predicted PTB. CONCLUSION: These findings represent an important step toward identifying key factors, which can be assessed rapidly before/after conception, to predict risk of PTB, and perhaps other adverse pregnancy outcomes. Quantifying these factors, before, or early in pregnancy, could identify women at risk of delivering preterm, pinpoint mechanisms/targets for intervention, and facilitate the development of interventions to prevent PTB. KEY POINTS: · Newly designed questionnaire used for rapid quantification of stress and psychosocial factors early during pregnancy.. · Deleterious factors predict increased preterm birth (PTB) risk.. · Protective factors predict decreased PTB risk..

In-silico generation of high-dimensional immune response data in patients using a deep neural network.

Technologies for single-cell profiling of the immune system have enabled researchers to extract rich interconnected networks of cellular abundance, phenotypical and functional cellular parameters. These studies can power machine learning approaches to understand the role of the immune system in various diseases. However, the performance of these approaches and the generalizability of the findings have been hindered by limited cohort sizes in translational studies, partially due to logistical demands and costs associated with longitudinal data collection in sufficiently large patient cohorts. An evolving challenge is the requirement for ever-increasing cohort sizes as the dimensionality of datasets grows. We propose a deep learning model derived from a novel pipeline of optimal temporal cell matching and overcomplete autoencoders that uses data from a small subset of patients to learn to forecast an entire patient's immune response in a high dimensional space from one timepoint to another. In our analysis of 1.08 million cells from patients pre- and post-surgical intervention, we demonstrate that the generated patient-specific data are qualitatively and quantitatively similar to real patient data by demonstrating fidelity, diversity, and usefulness.

Comparative Life Cycle Assessment of Chemical and Biocatalytic 2'3'-Cyclic GMP-AMP Synthesis.

Life cycle assessments (LCAs) can provide insights into the environmental impact of production processes. In this study, a comparative LCA was performed for the synthesis of 2'3'-cyclic GMP-AMP (2'3'-cGAMP) in an early development stage. The cyclic dinucleotide (CDN) is of interest for pharmaceutical applications such as cancer immunotherapy. CDNs can be synthesized either by enzymes or chemical catalysis. It is not known which of the routes is more sustainable as both routes have their advantages and disadvantages, such as a poor yield for the chemical synthesis and low titers for the biocatalytic synthesis. The synthesis routes were compared for the production of 200 g 2'3'-cGAMP based on laboratory data to assess the environmental impacts. The biocatalytic synthesis turned out to be superior to the chemical synthesis in all considered categories by at least one magnitude, for example, a global warming potential of 3055.6 kg CO2 equiv. for the enzymatic route and 56454.0 kg CO2 equiv. for the chemical synthesis, which is 18 times higher. This study demonstrates the value of assessment at an early development stage, when the choice between different routes is still possible.

Early prediction and longitudinal modeling of preeclampsia from multiomics.

Preeclampsia is a complex disease of pregnancy whose physiopathology remains unclear. We developed machine-learning models for early prediction of preeclampsia (first 16 weeks of pregnancy) and over gestation by analyzing six omics datasets from a longitudinal cohort of pregnant women. For early pregnancy, a prediction model using nine urine metabolites had the highest accuracy and was validated on an independent cohort (area under the receiver-operating characteristic curve [AUC] = 0.88, 95% confidence interval [CI] [0.76, 0.99] cross-validated; AUC = 0.83, 95% CI [0.62,1] validated). Univariate analysis demonstrated statistical significance of identified metabolites. An integrated multiomics model further improved accuracy (AUC = 0.94). Several biological pathways were identified including tryptophan, caffeine, and arachidonic acid metabolisms. Integration with immune cytometry data suggested novel associations between immune and proteomic dynamics. While further validation in a larger population is necessary, these encouraging results can serve as a basis for a simple, early diagnostic test for preeclampsia.